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Abstract 

The quantum dynamics of quarks and " gluons" and the scalar degrees of freedom associated with 
the non-linear regime of the non-linear gauge is derived. We discuss the subtleties in quantizing 
in a stochastic background. Then we show in detail that {SyMit^, 4''^', /")) j only depends on t^, 
thus effectively proving that the scalars (j)"" are non-propagating. Integrating out the scalars from 
(Sfermion) f leads to fermion wave function that declines exponentially. Finally, we derive the 
effective action of the "gluons" and fermions resulting from stochastic averaging. We show that it 
leads to a confining four-fermi interaction. 
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I. INTRODUCTION 



The non-linear gauge condition^ 

{d ■ D"^) {d ■ A'') = {D"^ ■d){d- A^) = {3^5"'' - ge^^'A^d^) (9 ■ A^) = 0, (1) 

and its physical consequences has been discussed by this author in a number of papers. This 
gauge condition is a natural generalization in Yang-Mills theory of the Coulomb gauge in 
Abelian theory. The reasons for this claim are: 

(a) In the Abelian limit, {d-D) and since is positive definite, the gauge-condition 
yields the Coulomb gauge. 

(b) The gauge condition has two regimes, the linear regime given by 9 ■ A" = 0, the 
Coulomb gauge, and the non-linear regime defined by 9 ■ A" = /"^ 7^ 0, which is also the 
zero mode of d ■ D. Thus the non-linear regime corresponds to the "Gribov" horizon on the 
surface d ■ A"- = /"(x). 

(c) The linear and non-linear regimes do not mix in the sense that field configurations in 
the non-linear regime cannot be gauge transformed to the Coulomb gauge and vice versa.— 

(d) If we consider the running of the coupling constant, we find that for the short- 
distance regime where g (asymptotic freedom phase), the non-linear gauge reduces to 
d'^{d- A"-) = 0, yielding the Coulomb gauge because d'^ is positive definite. This is consistent 
with the fact that transverse gluons are physical degrees of freedom at short-distances. 
However, as we increase the distance scale, g increases and before the running coupling 
becomes too large, where perturbation theory loses validity, the full non-linear character of 
the gauge condition becomes important where the relevant degrees of freedom are the scalars 
/"(x) and the "gluons" t^(a;), which decompose the Yang-Mills potential^ as 

A-^ = ^ ((5"^ + e^'j^ + rf){-d,f + tl). (2) 
(f ~l~ / ' /) ^ 

This means that the non-linear gauge continuously interpolates short-distance and large 

distance physics with their corresponding degrees of freedom. Since we consider quantum 

effects in the running of the coupling, the non-linear gauge is a quantum gauge condition. 

(e) Lastly, geometrically, configuration space analysis showed that non-linear gauge con- 
dition is a modification of the global orthogonal gauge condition, which does not exist in 
non- Abelian theories and gives the Coulomb gauge in Abelian theory.— 
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As for the physical consequences of the non-hnear regime of the non-hnear gauge, the 
following results had been established by the author: 

(a) The Yang-Mills action is quartic in and infinitely non-linear in The pure /" 
action has a factor and its kinetic term goes like [df^Y- Clearly, these hint of non- 
perturbative physics.^ 

(b) The pure dynamic shows that all spherically symmetric /"^(x) with x = {x^x^Y^"^ 
are classical configurations with zero field strength. The pure dynamics has a very broad 
minimum with zero action. Because of the infinite degeneracy in classical configurations, the 
author proposed to treat /" as a stochastic variable with a white-noise distribution. This 
resulted in the area law behaviour of the Wilson loop, which means a linear potential for 
static sources.'^ 

(c) If the stochastic treatment of the classical configuration /"(x) yields a linear potential 
between static sources, full quantum dynamics of shows equivalence to an 0(1,3) non- 
linear a model in 2D.® The proof made use of the Parisi-Sourlas mechanism.— Since the 
a model is 0(1,3) and not 0(4) resulting in a kinetic term with a wrong sign, the proof 
of confinement is purely formal. Furthermore, the proof of confinement only involves the 
scalars and not the "gluons" and quarks. 

(d) When we consider the classical dynamics of the "gluons" in the spherically symmetric 
background /"(x), it was shown that stochastically averaging /''(x) yielded a mass for the 
tf, and the loss of its self-interactions.— Thus, we have shown the mechanism for the mass 
gap. As for the loss of self-interactions, the result was also arrived at by Kondo.-^ 

(e) Finally, the author proposed the concept of limited gauge invariance and showed that 
in both the linear and non-linear regimes of the non-linear gauge we can define potentials 
that are gauge-invariant within their limited context.— 

In this paper, we will address the shortcomings noted in (c) by considering quarks and 
"gluons" in a stochastic background. We will decompose the scalar /"(x) via 

r(x) = r(x) + r(x). (3) 

We will note the subtleties in treating quantum and stochastic fluctuations. We will derive 
the effective action of quarks and "gluons" by averaging over the stochastic background. 
Finally, we will derive a confining non-local four-fermi interaction. 
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II. THE ACTION IN THE NON-LINEAR GAUGE 



We will consider SU{2) theory with the following action 

S = SyM + Sferrmon = J d^x{-^F;,F;, + ^l^l^ ^ij] . (4) 

Introducing the Fadeev-Popov trick of resolving unity using the gauge condition given by 
equation (0), we find the vacuum to vacuum functional 

1^(0) = j {dA1){di)){di))dete5{d ■ D{d ■ A))e-^ , (5) 

where 

gad ^ p ab g^^g . j^bd^ _ g^abc^df,{d ■ A^))Df. (6) 

Note that in the linear regime, where the coupling g is very weak and thus the smallest 
eigenvalue of the positive definite d"^ (in i?^) is not lowered to zero, {d ■ D) is still positive 
definite and the path-integral given by equation (jSJ reduces to the transverse gauge path- 
integral. 

As we further increase the coupling constant, we get to the Gribov horizon oi d ■ A"^ = Q 
surface. We are still on the Coulomb surface but the Fadeev-Popov operator is now singular. 
This regime had been discussed extensively in the late 1970's and 1980s. The conjecture is 
that restricting the transverse gauge path-integral to within the central Gribov regioni^ or 
the fundamental modular region results in confinementii. This author has a different view: 
confinement happens in the next stage. 

Further increasing the distance scale, we get to the regime where transverse gluons should 
no longer be relevant degrees. The reasons for this are the existence of the mass gap (massive 
vector fields are not transverse) and the fact that it is essentially impossible to get a confining 
interaction from the exchange of transverse gluons. Thus, we should be getting off the 
Coulomb surface and the non-linear gauge seems to be a natural gauge to consider. Here, 
we get to the Gribov horizon of the 9 ■ A" = /° 7^ surface where the zero mode of {d ■ D) 
is d ■ A°- itself. However, even though d ■ D is singular, the corresponding operator 9 given 
by equation Q is non-singulai^. The path-integral is still well-defined, presenting no need 
to deal with zero modes. 

In this non-linear regime, we can decompose in terms of the scalars and "gluons" 
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t° as given in equation (j2I). These new degrees of freedom satisfy the constraints 

a-r-^r = o (7) 

a ^ r abc , ^abd rd rc ^acd rd rb , ra j'd^dbc 

" (l + /-/> 

-r (1 + / ■ ns"^ - ni + / ■ f)5'^']djx = 0, (8) 

giving the same number of degrees of freedom as the original yl^(a;). Substituting equation 
(121) in equation (jH), we get 

S = [ d'x{{\z\f) + -Z{f) ■ LU\ t) + [2ZU) ■ Qif, t) + L{f, t) ■ L{f, t)] 
J g g 

+2gLif, t) ■ Qif, t) + g'Qif, t) ■ Qif, t)) + {i>t^,d,^P 

-zgi^j,Ty[R'^\f){-^d,f + <)])}, (9) 

where 

z;m') = x^'v)d,f%r, (10) 

L%{f,t) = R^\f){d,tl - d^) + Y^'%dJ% - djX), (11) 

Q;M,t) = T'^''VKt:, (12) 

X'^^^f) = + 2/ ■ ne'^^^ + 25'^V - 25"7^ 

_^^^abdjdjc _ ^^acdjdjb ^ (:bcdfajd^^ ^^3^ 

= ^^^(5"^ + e^^T + rf), (14) 

(1 + /-/) 

Y^"'^ = - — \^A-{f- + (1 + / ■ f)n'' - (1 - / ■ fW'^f 

(1 + / ■ fr 

j^^^cadjdjb _ 2pff^ + ^^bdjdjc ^ jagbcdjdj^ ^^5^ 

j.abc ^ + (1 + / ■ f)f5^^ - (1 + / ■ f)r5^' 

(1 + / ■ fr 

I ^abd jd jc _j_ ja^cd jd _j_ ^acd jd y6j (16) 

The pure f"' dynamics given by Z^ has a class of classical configurations with zero field 
strength: the spherically symmetric /'^(x), with x = {x^XyY^"^ (we are in R^). This follows 
from the anti-symmetry of X"-^'^ with respect to the last two indices and that for spherically 
symmetric /"(x), 9^/°(x) = Since we will be expanding about a class of classical 

configurations, the author proposed to treat the spherically symmetric f°'{x) stochastically 
with a white-noise distribution. 
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What will be done in the remainder of this section is to substitute the background de- 
composition given by equation Q in the action given by equations and (fTUj) to (fTBjl . 
The resulting action is quantic in and infinitely non-linear in (f)"" with coefficients that are 
ftmctions of f°'{x). The result is (see Appendix A). 

oo 
n=3 

^ X X ^ ^ 

+ t) ■ g(/, t) + (7'g(/, t)g(/, t)]^^^^' (a^i) ■ ■ ■ 0"" 

°° 1 pab/' f 1 1 

+ g ;^ ,/.(.0 . . .//.■(.„) '-7^ + + '^>^°''--'' ■ ■ ■ (17) 

The first term bracketted with [ ] comes from Z'^{f) where we made use of Z^^{f) = 0. 
The second bracketted terms have been discussed in a previous paper^, where the author 
presented the classical dynamics of "gluons" in a stochastic vacuum. The third shows the 
interactions between the "gluons" and scalars 0". The fourth gives the fermion action, with 
the interaction between the fermions and "gluons" and scalars in a classical background 
f°'{x). Equation (|T7j) looks very complicated where the scalars are infinitely non-linear and 
the background f"'{x) is rather involved. However, this complicated action is dramatically 
simplified upon stochastic averaging. 

Introduce the white noise distribution 

P[~f] = Afexp.i-j ns)ns)ds} (18) 

where ^ is the scale when non-perturbative physics is important. From the running of the 
coupling, we must have £ ~ ^qcd- When we average equation (fTTj) using equation (fTH|). a 
tremendous simphfication results (see Appendix B) 

{s)f = J{dr{x))s{t;,r,^,tj~)p[f] 
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From equation (fTT^ . we see that the scalars 0" do not have a kinetic term; they do not 
propagate. 

When we get the equation of motion for 0", we find 

+ A-y^i^ = 0- (20) 

The spherically symmetric solution of equation ()2()|1 is 

tpix) = iJoexp.{--^g^x}, (21) 

which is an exponentially vanishing fermion field with effective length = ~ AqJ^j^. This 
behaviour is a signal of confinement. 



III. QUANTUM FIELD THEORY IN A STOCHASTIC BACKGROUND 

We would like to quantize t^,?/','^ and 0" in the presence of the stochastic background 
In essence, we need to deal with quantum fluctuations on top of classical stochastic 
variables. This poses an ambiguity due to the ordering of how these two fluctuations are 
taken into account. One sequence is to consider quantum field theory in a stochastic back- 
ground and then do the stochastic averaging. The other sequence is first to do the stochastic 
averaging on the classical action and then quantize the resulting theory. 

Let us collectively represent the fields t^, tp, ip and (j)'^ by the field $. The stochastic classi- 
cal background is still given as /. We write the action of $ in a background /, which includes 
gauge-fixing and the "Fadeev- Popov" determinant as 5'[$,/] = / (i''a;£($, /). Define 

= ^e//[$]. (22) 

The complete generating functional from which we compute n-point functions is 

Weff[J] = [ (t/$)e-^^^^[*l-^^'"''* (23) 
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From this generating functional, we compute the n-point function by 

G„(xi, ■ ■ ■ , = ^j^^^yHD^^^y (24) 

On the other hand, we can first consider the complete generating functional in the pres- 
ence of the background /. This is 

W[J- f] = j (rf$)e-^[*'^'l--^'^'"-'*. (25) 

Doing stochastic averaging, we find 

WlffiJ] = {W[JJ])f 

(rf<l')(e-^[*'^'l)/e~^'''"-^* (26) 

Expanding the exponential and taking the stochastic average of each term, we find 

(e-^[*Jl)^- = e-^^//[*l, (27) 

where 

+ ■■■ (28) 

The terms £($, /; x)C{^, /; y), etc., are correlated points, which arise because the deriva- 
^ ^ ' 

five of the white noise f{x) is " smoothened-out" via 

^ - I ^^^^ 

n 

Using equations (P7j) and in (|^. we find that W^jrj[J] ^ VFe//[J]. The n-point function 
before the stochastic averaging is given by 

^"^"^•••'"-^^ = ^J(a:0---Mx„) - 



Obviously, 



(C(x x-f)) -( ^"^[-^'^'J ) 



6J{xi) ■ ■ ■ 6J{Xn) 

G;(xi,---,a;„) (31) 
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However, it is clear that the n-point function that appear in equation ()31|) is not equal to 
the n-point function given in equation (j211). 

Thus, we see that the method of first stochastically averaging the classical action in the 
presence of the background / and then quantizing the theory is not the same as one that 
first quantizes theory in the presence of / (see equation ((201)) then does the stochastic 
averaging. 

There is another subtlety in quantizing the theory even within the framework of doing 
the stochastic averaging only after quantizing the theory in the presence of /. Let us begin 
with equation ()25|) . which represents the complete generating functional in the presence of 
the stochastic background /. This is diagrammatically represented by 

W[JJ]^@ (32) 

while the n-point function given by equation (jHDI) is represented by 



Gn{Xi ■■■ ,Xn]f) = <^ (33) 

Equation (31) says that 



(34) 

where the RHS of this diagrammatic expression represents the full n-point function from 

Following the usual field theory prescription, we define the connected Green function 
generating functional via 

W\J- f] = e'^t-^'^'l (35) 
Diagrammatically, this is represented by 

(D = Z[J-J]. (36) 

We compute the connected n-point Green function via 

G^ix, ■■■ x-f)- (37) 
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It is clear from equation ()34|) that 



{Z[J;f])^-=-^{inW[J;f])^- 

^^inW^^f[J]=Z'^ff[J] (38) 

where the last hne of equation (|HH|) follows from equation ^2^. From equation (|H7j). it follows 
that 

/))/7^G,}/,n(^i,---,a;„), (39) 
where the RHS is evaluated using equation (f^ . Diagrammatically, this is represented by 



7^ 



(40) 

If this non-equivalence is true for the connected n-point function it is easy to show that 
it is also true for the one-particle- irreducible functions derived from the r[$,/] and r[$], 
which are derived via Legendre transformation. Define the "classical" fields 

I. -^SiAH (AO) 

where equation makes use of equations (^Hjl and and it is obvious that 

The effective action, which is also the one-particle-irreducible generating functional is 
defined by 

r[|., /] = z[j, f] - J (fx^CfU (43) 

Kffi^eff] = Z'^ffiJ] - I d'xJ^^ff. (44) 

where it is implied that equations (PT|) and (jl^ are inverted to solve for the sources in terms 
of the "classical" fields. Again, it is clear that 

(f[|.,/])^-^f^^/<|.e//]. (45) 
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Here, we see that although there is ambiguity when stochastic averaging is done at the 
level of connected and one-particle-irreducible generating functionals and n-point functions, 
there is no such ambiguity at the level of the full generating functional and full n-point 
functions. Since the S-matrix is expressed in terms of the full n-point functions, the quantum 
theory in a stochastic background is consistent and well-defined as long as we make use of 
the action given by equation (j^ . 



IV. THE EFFECTIVE QUANTUM DYNAMICS 

We begin with the path-integral in the non-linear gauge given by 

w = j{dti){dn{dmmd-t''-^,nKp'') 

■ det~'^{l + f- f)det9eXp.{-{SYM + Sfermion)}- (46) 

Let us rewrite the delta functionals by 

■ - = det[ ■ r - -^n) (47) 

gP' (1 + /-/)^ (1 + /-/)^ 9^' 

dip") = detl i^^]5( ^^^p1 (48) 

(1 + /-/)^ {1 + f-fr 

where j and k are positive integers. The reason for equations ()47|) and ()48|) will become clear 
later. Expressing the determinants in terms of ghosts, the path-integral can be written as 

W = j {dt'l){dP){di)){di)){du''){dvf')exp.{-S'}, (49) 



where 



S' — SyM + Sfermion + Sgf + Sghosts, (50) 

= [ d^x{{-) — l^^id ■ e - -^ry + 1 — L^(p«)2} (51) 



Snhosts = / d^xu"" „ e^'^u^ (52) 

Implementing the background decomposition given by equation (jS)) in the path-integral, 
we get the vacuum to vacuum functional in the background /"(x), i.e., 

W[P] = [ idtl){d^''){dtP){di'){du''du'')exp.{~S'}, (53) 
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where S'{r^,ilj,iljJ'' = + As equations ^ and ^ show, as far 

as the full generating functional and the n-point Greens functions are concerned, replacing 
S' by S'^jj- (see equation (jlHl)) in equation is equivalent to computing physical processes 
with and then doing the stochastic averages at the end. Now we will see the reason 

for equations (jTTjl and pHjl . Equation says we have to evaluate 

{S') f = {SyAi) + {Sf ermion) + {Sgf) + {Sghosts) (54) 

/d xd yii^CyM ~\~ ^fermion ~\~ ^gf ~\~ ^ ghosts) x{,^Y M ~\~ ^fermion ~\~ ^gf ~\~ ^ghosts)y) f ; (5^) 
V ' 

etc. 

The stochastic averages involve integrals of the form 

0, for m < n 

lim(7r~^/^o-+^/^) / -; — e"'^'' (ir = < non - zero, finite, form = n + l (56) 

diverges, form > 0, n + 2. 

Because of the - — jT^ term that goes with (9 ■ 1?/")^, when we expand using the 

decomposition given by equation Q and by suitable choice of j, the resulting integral will 
always yield zero. The same thing is true with the term ~ (p^) and the ghost term. Thus, 
i^af) f ~ {^ghosts) / = 0- -^'^^ {^') f yields just the term given by equation (fT^ . 

As for equation all the correlated terms involving Cgf and Cghosts vanish for the 

same reason as above. This means that 5*^^^ does not involve any ghosts m", u"" and we can 
just lump the ghosts measure with the normalization of the path-integral. 

The fermion-fermion correlated term will yield the following non-local, four fermi inter- 
action 

2 r 

NLFF = ^J d'xd'y{i,^,T'^ij)Ml{^)At{y))f{i^luf'^'i^)y (57) 
We will evaluate the stochastic average by making use of 



From equation we must have 



&',dy{AU^)Aliy))f = ^5(x -y) + ..., (59) 
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where x = {x^x^y/'^ and y = (y^y^)^/^. Equation imphes that 

it for a sph 

and -^l^ — y\ = 5{x — y). Substituting equation (jUOj) in NLFF, we find 



{Al{x)K{y)) = -y\ + ... (60) 

The equivalence follows from the fact that for a spherically symmetric function, = -^-^ 



NLFF = ^(i) j d^xd%{^r^,^,T^ijUx - y|(^7.^,T», + ... (61) 

where 

= {sin9isin62sin(j), sin9isin92C0S(f), sin9iCos92, cos9i), (62) 

i.e., r]fj_ represents the unit vectors in 4D spherical coordinates. If the fermion field is spher- 
ically symmetric, the angular integration does not vanish only when r7^(x) = ±ri^{y), i.e., 
the 4D vectors are collinear. Using 

dVL^'q^T]^ = —5^^ (63) 

and J dVti^ = J^^ sin^ 9i sm92d9id92d(f) = 27r^, we find that we can write the equation 
dnH) as 

NLFF = i(i) J d'xdS{i'l,T'^i^US-mi,T^i^)y + " " " (64) 

where the points x and y in i?^ must be co-linear, suggesting "fiux-tube" configuration. This 
is surprising because the mechanism involves spherically symmetric, stochastic f°'{x) yet the 
four-fermi interaction yields a co-linear configuration with linear interaction. Furthermore, 
if we neglect the extra terms in equation ()64|1 (which will involve quadratic terms in cj)""), 
the 0" integral will yield the following delta functional 

5(9,(V^7mT» - ^^(^7,T'^^)). (65) 

As already discussed in section II, this implies equation ()20j) yielding fermion solutions given 
by equation (PT|) . which points to fermions with effective length ~ ^^cd ^"^^ with a linear 
potential as given by equation (j64p . This definitely shows quark confinement and that when 
we try to pull a fermion from inside a hadron, a colorless quark-anti quark pair (a meson), 
will be formed because of the linear potential. This is the picture of "string" breaking, 
although a string was never invoked in our mechanism (spherically symmetric, stochastic 

~nx)). 
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In summary, the effective theory for quarks and "gluons" under the approximations we 
make, is given by the path-integral 

X ea:p.{-%/(t^,V,^)} (66) 

where 

+ i(i) j d\d'y{i,^,TyUx-y\{i,^,Ty)y. (67) 

V. CONCLUSION 

We have derived an effective dynamics for quarks and "gluons" as given in equations (jMj) 
and (jHTj) . The effective action clearly shows a mass gap, i.e., the "gluons" acquired a mass 
and the quarks are confined. The mechanism for all these is the spherically symmetric /"(x) 
treated as a stochastic variable. These vacuum configurations arise from the non-linear 
regime of the non-linear gauge, which we claim is the natural generalization in Yang-Mills 
theory of the Coulomb gauge in Abelian theory. 
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APPENDIX A 



Here, we will derive equation (fTTj). The starting point is equation Q. Using equation 
in equations (fTUI) to (fTB|) . we get 



X (9,/'' + 9^0'')(9.r+9.0^) (Ai) 

Since X^^^{f)dJ^dJ'' = 0, we find that 

K^lif) = 2X'^'^{f)X'^^\f)dJ'dJ'6,, - 2X^''^{f)X^^\f)dJ%r (A2) 

rii^cde _ XX^fd _____ 

M""!/) = -jj^ifV-jjjrimWdJ^dJ^ (A3) 

AX'^^'/' _ _ _ 

iVf (/) = AX'^^^dJ'dJ%ff (A4) 

The rest of the terms in the first bracket of equation p7|l can be arrived at by considering 
the second and higher functional derivatives of X"''^ 

As for the second bracketted terms in equation ()17p. these had been derived in reference 
[8]. The third bracketted terms, which give the scalars-"gluons" interaction in a background 
makes use of equations to (fTH|l . This is an infinite series of complicated terms. The 
important point is that each term has sufficient powers of (1 + / ■ /) in the denominator, 
which makes it vanish upon stochastic averaging. 

As for the fourth bracketted terms, these arise from 

= R-''{f + ^){^dJ''+-d,<f>' + tl) 
9 9 

{Uj' + U,<p' + tl). (A5) 
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APPENDIX B 



We will derive equation (fT^ from equation (fT7|). We will make use of equation (fTHj) and 
(jSni). Let us evaluate in detail one stochastic average. 

= AfJ {dnx^''\f)X^^\f){dJ''dJ^5,,-dJ''dJ^) 
X exp.{-j ds~r{s)r{s)}. (Bl) 
Using equation (|^ and equation (fT^. we find that we need 

(/'(^ + -))/ = (B2) 



+ = i*''ito(-) (B3) 

where a = The stochastic average at point x, which goes with ^ is of the form 

+3e'='^-^a;-^x" - Se^^^x-^x'^ + e'^'^-^x^x^] [-(1 + 2r'^)e''^^ 
+25''x^ - 2(5"''x" + ^e'^'x^x^ - ^e'^'^x^x' + e'^^'^x^x"] 

xe-"'^ (B4) 



Since the numerator only has 6 powers of r at most, while the denominator has 8, the 
integral is at best 



^-3/2^+3/2 



/•oo -crr^ 

/ ^ dr = 7r-^/V-3/2[l - $(a3/2)-e'^] (B5) 

Jq l + r^ 2 



where $ is the error function (not to be confused with the field that collectively represents 
0**, 0, ip in section III) with expansion 



Thus, this integral behaves, at best, like cr"*"^/^, as a — > 0. Combining this with the ^ factor, 
we find this term to vanish like a^/^. There is another term in {K^^) j with /'^(x)/^(x) of 
the derivatives, which will combine with the rest given in (B.4). But this term will yield a 
stochastic average which will vanish like a as a — > 0. Thus, we find 

(^S(/))/ = 0. (B7) 
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Using similar analysis, it is trivial to show that all the other stochastic averages in the first 
and 3rd brackets of equation ()17|) vanish. The only non-vanishing stochastic averages are 

{R''\f))f- = Is^' (B8) 
{Y^')j = 6^\-){-) (B9) 

These averages yield the first and second terms of equation (fT^ . 
Finally, from equation (A. 5), we get 

iK)f-t; + Ur-l(jf-^r. (BIO) 
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